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Abstract 

We have combined high resolution VLBI Space Observatory Programme (VSOP) data at 1.6 and 4.8 
GHz with Very Long Baseline Array (VLBA) data at higher frequencies and with similar resolutions 
to study the spectral properties of the core of M87 with milliarcsecond resolution. The VSOP data 
allow a more accurate measurement of the turn-over frequency, and hence more reliable determination of 
associated physical parameters of the source. Comparison of the images with previously published images 
yields no evidence for significant motion of components in the parsec-scale jet. In addition, the brightness 
temperatures obtained from model-fits to the core are well below the inverse Compton limit, suggesting 
the radio emission we are observing is not strongly Doppler boosted. 
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1. Introduction 

The identification of the radio source Virgo A with 
the giant elliptical galaxy M87 was one of the first such 
associations to be made and demonstrated the gains in 
angular resolution obtainable from interferometric obser- 
vations (Bolton et al. 1949). M87 (NGC4486, 3C274, 
J1230-I-1223) lies near the center of the Virgo Cluster and 
was already well known due to its "curious straight ray" 
(Curtis 1918). This jet extends over 20" at optical (e.g.. 
Sparks et al. 1996) and X-ray (e.g., Marshall et al. 2002) 
energies and has been well studied at radio wavelengths 
(e.g., Owen et al. 1989). VLA observations at 90 cm have 
revealed the jet extends considerably further, appearing 
to flow into two large bubbles in the intergalactic medium 
(Owen et al. 2000). 

Hubble Space Telescope (HST) spectroscopic observa- 
tions indicate the presence of supermassive black hole 
in the core of the galaxy, with 3 x 10^ Mq contained 
within the central 3pc (Harms et al. 1994) and the 
process of accretion of material towards the core is be- 
lieved to power the jet. At a distance of only 16.75 Mpc 
(Whitmore et al. 1995), high angular resolution stud- 
ies correspond to the highest linear resolutions obtain- 
able for AGN, with 1 mas corresponding to 0.08 pc. The 
inner parsec-scale jet has been studied at increasingly 
higher radio frequencies, revealing that jet coUimation be- 
gins within several tens of Schwarzschild radii of the core 
but is not complete until several thousand Schwarzschild 
radii (Junor et al. 1999; Ly et al. 2004a; Krichbaum 
et al. 2004; Krichbaum et al. 2005). A comparison of 
Chandra and HST data indicates that the X-ray flux dom- 
inates closer to the core but decreases relative to the 
optical flux with increasing core distance (Marshall et 
al. 2002). A detection of M87 has also been reported at 
TeV energies from a region which includes the core and 
inner jet, although systematic pointing errors prevented 



the identification of the gamma-ray production site with 
any feature at other wavelengths (Aharonian et al. 2003). 
The TeV detection has, however, yet to be confirmed (Le 
Bohec et al. 2004). 

The study of motions of features in the jet has yielded 
a number of surprises. Early VLBI observations revealed 
that apparently superluminal motions (i.e., Wapp > c) 
were common in radio-loud active galactic nuclei (AGN), 
however the first multi-epoch 1.6 GHz VLBI observa- 
tions of M87 suggested that the parsec-scale jet displayed 
sub-luminal (i.e., Wapp < c), with an apparent speed of 
l.l±0.3masyr~^, corresponding to 0.28±0.08c (Rcid ct 
al. 1989). Observations with the VLA at 15 GHz at six 
epochs spanning 9 years resulted in the detection of ap- 
parent motions between 0.5 c and ~ c for the bright knots 
over 1 kpc from the core, but apparent speeds of ~ 2c for 
two features ^ 200 pc from the core (Biretta et al. 1995). 
Further VLBI observations at 22 GHz yielded only up- 
per limits, with apparent speeds of ^ 0.3c for components 
within 4 mas of the core (Junor & Biretta 1995) and 
the analysis of further 1.6 GHz VLBI images suggested 
features out to 160 mas were stationary (—0.03 ± 0.02c, 
Biretta & Junor 1995). 

In contrast, a series of five HST images resulted in 
the measurement of apparent speeds of 2.6c-6c for com- 
ponents between 0.87" and 6" of the core (Biretta et 
al. 1999). An apparent speed of 0.63±0.27c was mea- 
sured for the innermost resolvable component, 160 mas 
from the core (Biretta et al. 1999). More recently, 
Kellermann et al. (2004) have reported an apparent speed 
of 0.14±0.07masyr~\ or 0.04±0.02c for a component 
6 mas from the core, based on 9 epochs of 15 GHz observa- 
tions over the six years from 1995 to 2001. A component 
detected between 1 and 2 mas from the core showed no 
net motion over the same time range. Finally, at 43 GHz 
Ly ct al. (2004b) have reported lower limits on the appar- 
ent speed of 0.25 c to 0.40 c for the jet, and 0.17 c for the 
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counter-jet. 

We report here on contemporaneous observations of 
M87 made as part of the VLBI Space Observatory 
Programme (VSOP; Hirabayashi et al. 1998; Hirabayashi 
et al. 2000) and the Very Long Baselme Array (VLB A). 
The multi- frequency data allows the evolution of the spec- 
tral index with core distance to be studied, and compari- 
son with previously published images enables component 
motions to be examined. 

2. The Observations 

A series of VSOP observations of M87 have been made 
at 1.6 GHz and 4.8 GHz over the lifetime of the mission 
(see, e.g., Junor et al. 2000; Biretta et al. 2002). We have 
selected from the VSOP data archive 1.6 and 4.8 GHz ob- 
servations made over a three day period in March 2000, 
and combined them with VLBA data from observations 
made three weeks later, to minimize the effects of source 
variability. Details of the observations presented here are 
given in table 1. 

The VSOP observations were made in the standard 
VSOP observing mode, with two 16 MHz bandwidth 
channels of two-bit sampled, left circular polarization 
(LCP) data recorded at all stations (Hirabayashi et 
al. 2000). HALCA data were received and recorded at the 
Tidbinbilla, Goldstone, Robledo and Green Bank tracking 
stations. The ten telescopes of the Very Long Baseline 
Array (VLBA) co-observed at both 1.6 and 4.8 GHz, and 
in addition the Tidbinbilla and Robledo 70 m telescopes 
participated in the 1.6 GHz observation. 

The VLBA observations on 2000 April 7 recorded 8 IFs 
(both polarisations) of 8 MHz bandwidth alternating be- 
tween 15 and 22 GHz approximately every twenty minutes. 
North Liberty was found to have poor signal to noise and 
was excluded. The 43 GHz observation on 2000 April 9 
was made with 8 IFs of 8 MHz bandwidth, and included 
the phased VLA, but not Hancock. 

3. Data Reduction 

The data were reduced in AIPS^ using the standard 
VSOP and VLBA methods, and the VLBA pipeline 
scripts where possible. The VLBA antenna system tem- 
peratures are recorded within the data files. The nomi- 
nal system temperatures for the respective observing fre- 
quencies were used for HALCA (Hirabayashi et al. 2000) 
and also DSN telescopes, and the VLA antenna tempera- 
tures were scaled to 3C286 (also observed in the exper- 
iment). All of these were applied to the data via the 
task ANTAB. After this initial amplitude calibration, de- 
lay and phase calibration was undertaken with fring, be- 
fore exporting the data for deconvolution and imaging in 
Difmap (Shepherd 1997). For further image analysis we 
used MIRIAD (Sault et al. 1995) for alignment, smooth- 
ing and spectral index calculation. Finally we used Octave 

^ Astronomical Image Processing Software developed and main- 
tained by the National Radio Astronomy Observatory. 



(Eaton 2002 ) to calculate the turn-over frequency using 
the minimisation routines therein. 

4. Images 

We present images of M87 at 1.6, 4.8, 15, 22 and 
43 GHz at their maximum resolutions in figure 1. The 
full-resolution VSOP image at 1.6 GHz has a resolution 
of 1.23x0. 96 mas at a position angle (PA) of 77°. It was 
made with uniform weight binning and no amplitude error 
weighting to achieve maximum resolution at the expense 
of signal to noise. The peak brightness is 0.20Jy/beam 
and the RMS is l.lmJy/beam. The low-level emission 
surrounding the core perpendicular to the jet axis in the 
1.6 GHz image has not been reported previously. It is not 
seen in higher frequency images, consistent with this weak 
extended emission having a steep spectral index. This 
emission would not have been visible in lower resolution 
images such as those in Rcid et al. (1989) and Reid (1998). 

The 4.8 GHz VSOP image has a resolution of 0.7 x 0.5 
mas at a PA of —82°. It was made with uniform weight- 
ing and an amplitude error weighting of —1. The peak 
brightness is 0.14 Jy/beam with an RMS of 0.2mJy/beam. 
The image is similar to the VSOP images in Junor et 
al. (2000) and Biretta et al. (2002). For analysis of spec- 
tral variations it is necessary to convolve all images with 
the same beam-size: At 4.8 GHz the smoothed image has 
a resolution which matches the 1.6 GHz image, a peak of 
0.14 Jy/beam and an RMS of 0.3mJy/beam. 

The 15 GHz VLBA image has a resolution of 0.9 x 0.5 
mas at a PA of 2°. It was made with uniform weight- 
ing and an amplitude error weighting of —1. The peak 
brightness is 0.56 Jy/beam with an RMS of 2mJy/beam. 
The smoothed image has a resolution which matches the 
1.6 GHz image, a peak of 0.83 Jy/beam and a RMS of 
2mJy/beam. The image is similar to those made as part 
of the VLBA 2 cm survey (Kellermann et al. 2004). 

The 22 GHz VLBA image has a resolution of 0.6 x 0.3 
mas at a PA of —1°. It was made with uniform weight- 
ing and an amplitude error weighting of —1. The peak 
brightness is 0.44 Jy/beam with an RMS of 1.5mJy/beam. 
The smoothed image has a resolution which matches the 
1.6 GHz image, a peak of 0.86 Jy/beam and a RMS of 
3mJy/beam. This can be compared with the 22 GHz im- 
age of Junor & Biretta (1995). The weak feature about 
1 mas to the North-East of the core could belong to the 
counter-jet and matches the direction (but not distance) 
of a similar component seen at 86 GHz (Krichbaum et 
al. 2004). 

Our 43 GHz image has a resolution of 0.3 x 0.2 mas 
at a PA of —4°. It was made with uniform weight- 
ing and an amplitude error weighting of —1. The peak 
brightness is 0.53 Jy/beam with an RMS of 1.5mJy/beam. 
The smoothed image has a resolution which matches the 
1.6 GHz image, a peak of 0.63 Jy/beam and a RMS of 
5 m Jy/beam. The peak brightness is significantly higher 
than that reported by Junor et al. (1999) for March 1999 
but less than that in Ly et al. (2004a) observed in October 
2001. The structure is consistent with those of Junor et 
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Table 1. Observing Sessions 



Observation code 


Frequency 


Date 


Ground Telescopes* 




[GHz] 






w022a7 


1.646 


2000 Mar 22 


HALCA, VLBA, Ti, Ro 


w040a5 


4.812 


2000 Mar 20 


HALCA, VLBA 


bjOSla 


15.365 


2000 Apr 07 


VLBA 


bjOSla 


22.233 


2000 Apr 07 


VLBA 


bjOSlb 


43.217 


2000 Apr 09 


VLA, VLBA-Hn 



* Ground telescopes: Hn = Hancock 25 m; Ti = Tidbinbilla 70 m; Ro = Robledo 70 m; 
VLA = phased VLA; VLBA = full VLBA. 



Table 2. Observation Parameters 





Resolution 


Residual 


Residual noise 


Core Brightness 


Frequency 


(beam size) 


noise level 


for 1.23x0.96 beam 


Temperature 


[GHz] 


[mas] 


[mjy/beam] 


[mJy/beam] 


[K/lQiO] 


1.646 


1.2 X 1.0 


1.1 


1.1 


2.8 


4.812 


0.7x0.5 


0.2 


0.3 


3.2 


15.365 


0.9 X 0.5 


2.0 


2 


1.8 


22.233 


0.6 X 0.3 


1.5 


3 


1.4 


43.217 


0.3 X 0.2 


1.5 


5 


5.7 



al. (1999); Biretta et al. (2002); Ly et al. (2004a), although 
with a lower signal to noise. Overlaid on our image (fig- 
ure le) are the opening angles identified by Junor et al. 
(1999). 

For each of these high resolution images we have calcu- 
lated the core brightness temperature in the source rest 
frame from the core integrated flux and deconvolved size 
and these are reported in table 2. We have taken the 
Eastern- most component as the core, which is also the 
brightest for all but the 1.6 GHz image. In the 1.6 GHz 
image the brightness temperature of the peak component 
is 11 X 10^° K. These values are all well below the in- 
verse Compton limit and therefore do not require Doppler 
boosting. 

5. Jet Component Motions 

As mentioned previously, M87 has been well studied 
using the VLBI technique and there are many published 
images with which we can compare those presented here. 

First, we have re-imaged the 1.6 GHz data with only the 
VLBA telescopes, using a Gaussian (u, v) taper falling to 
0.5 at 50 MA. The resulting image, shown in figure 2a, 
can be directly compared with the 1.6 GHz image of Reid 
et al. (1989) from epoch 1984.26, and also with those of 
Biretta & Junor (1995) at epochs 1988.41 and 1992.47. 
The similarities between the images are startling. Three 
jet features were identified by Reid et al. (1989): N2, at 
a core distance of ^20 mas, Nl at '^65 mas, and M at 
^^170 mas. Component M was subsequently labeled L in 
Biretta & Junor (1995) and Biretta et al. (1999). AU 
three components are visible in figure 2, and appear to 
have remained stationary over the 15.96 years spanned 



by these observations (with Nl and N2 having been first 
identified in lower sensitivity observations several years 
earlier). We place the positions of knots N2 and Nl at 
19.4 and 64.7mas along the jet in the images from Reid 
et al. (1989) (epoch 1984.26), to better than 2 mas. From 
our observations we place N2 and Nl at 18.4 and 62.9 mas 
in epoch 2000.22. We can therefore confirm the upper 
limits of -0.02 ± 0.03c and -0.03 ± 0.03c for N2 and Nl 
respectively. 

Biretta et al. (1999) reported a motion for knot L of 
2.48±1.06masyr~^ based on HST observations. This 
would suggest motion of at least 22.7 mas over the 15.96 
years, however the peak of the radio emission in knot L 
shows no evidence of having moved at all. In addition, 
the "gap" or minimum in the jet emission ^llOmas from 
the core does not appear to have moved. 

Closer to the core, we can compare our VSOP 4.8 GHz 
image with the VSOP image at epoch 1997.97 (Junor et 
al. 2000; Biretta et al. 2002). In addition to N2, enhance- 
ments in the jet emission can be clearly seen at ~13mas 
and ~6 mas, the latter corresponding to component B of 
Kellermann et al. (2004). These have displayed no dis- 
cernible motion over the intervening 2.25 years, although 
the speed reported for component B by Kellermann et 
al. (2004) of 0.14±0.07masyr"i cannot be ruled out. 

These observations strongly support the theory that 
these components in the radio jet mark the sites of stand- 
ing shocks. Filamentary features, limb-brightening, and 
side-to-side oscillations have been reported on a number of 
occasions (e.g., Owen et al. 1989; Reid et al. 1989; Biretta 
et al. 2002) which, as these authors have pointed out, is 
consistent with a fast underlying fiow. The visible radio 
emission arises from a less strongly beamed interaction be- 
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tween the inner jet and the surrounding sheath. The emis- 
sion from the sheath must still be Doppler boosted to some 
extent as otherwise the counter-jet would be expected to 
be easily detectable (e.g., Junor & Biretta 1995). It is 
possible to estimate the required Doppler factor from the 
ratio of the jet brightness to the noise level. Following 
the analysis of Jones et al. (1996) we use the following 
expression for smooth continuous jet sheath flow: 

/ l + /3cos^ y+" 

where R is the flux ratio between the jet and counter-jet 
(the 3cr upper limit in this case), /3 is the speed of the jet 
sheath in terms of the speed of light and 9 is the jet angle 
to the line of sight. The spectral index, a, has the opposite 
sign in this paper than in the reference, as discussed in 
section 6.1. At the peak of the 1.6 GHz emission the 
brightness is 0.2Jy/beam, a is 1.9 ±0.2 and the image 
noise is l.lmJy/beam. For 6 we use the value of 19tlo° 
of Biretta et al. (1999). We find R > 60, therefore the 
sheath has P > 0.5. 
The doppler factor, 

5=[(l-/32)-i/2(i_/3cos0)]-i 

implied by these values is therefore S > 1.6. 

The discrepancy between HST-based speed for knot L 
with the stationary radio features may perhaps be re- 
lated to the phenomena observed in the HST observations 
for the component HST-1 East, 870 mas from the core. 
This knot, which has a speed of 0.84 c, appears to launch 
new optical features which move at speeds of up to ~ 6c 
(Biretta et al. 1999). Possibly the stationary radio fea- 
ture 170 mas from the core is the site from which optical 
components like knot L are launched (with, in this case, 
sub- luminal speeds). We note that there are no published 
images of M87 at 1.6 GHz from the period of the HST ob- 
servations (1994.59 to 1998.55) to determine whether the 
HST feature has a radio counterpart. 

Biretta et al. (1999) suggest that the ~ 6c motions re- 
flect the underlying speed of the jet. It is clear that fea- 
tures in the radio jet between ~7 and 160 mas appear 
stationary, in contrast to the lower limits of 1 c to 1.6 c 
motions recently reported from observations at 43 GHz 
(Ly et al. 2004b) . It will be interesting to see whether the 
actual speeds close to core are comparable to the under- 
lying jet speed inferred from the HST observations. 

Finally, as described in more detail in the next section, 
we detect a component '^2 mas from the core in the 1.6, 
4.8, 15 and 22 GHz images. Precise measurement of the 
core distance is complicated by frequency dependent opac- 
ity effects and, potentially, by the time intervals between 
observations, if the lower limits to motions derived by Ly 
et al. (2004b) still apply at these core distances. We note 
that this core distance is similar to that of the compo- 
nent C of Kellermann et al. (2004). As noted earlier, com- 
ponent C was detected between 1 and 2 mas from the core 
by Kellermann et al., but showed no net motion over the 
six years. It has been suggested that "stationary" compo- 
nents may be standing recoUimation shocks, which could 



appear to move back and forth to some degree as other 
jet components move through them (sec, e.g., Jorstad et 
al. (2001) and references therein): it would be interest- 
ing to monitor the movements of this feature to deter- 
mine whether there is any relationship with the motions 
of the moving jet components observed closer to the core 
at 43 GHz. 

6. Image Alignment 

Self-calibrating VLBI data loses the absolute position 
information, and produces images centred on the strongest 
emission, so re-alignment is required when comparing im- 
ages at different frequencies. There are various approaches 
to re-aligning observations: we chose the simplest, which 
also produces essentially identical answers to more com- 
plex methods. 

The jet in M87 is straight, so does not allow the align- 
ment on a frequency-independent feature such as a bend. 
The region closest to the core could be approximately de- 
composed into two components between 1.6 and 22-GHz, 
however, so it was possible to identify and align on the core 
and innermost jet component. The latter, which is '^2 mas 
from the core, is the brightest feature at 1.6 GHz, with 
the core being brighter at 4.8 GHz and above. Because of 
the poor representation, and therefore positional accuracy, 
of the Gaussian components (particularly for the 1.6 GHz 
image) we abandoned unconstrained image alignment. We 
finalised on identifying 1-D Gaussian components along 
the jet axis. This allowed us to produce images aligned 
on the maximum of the core, i.e., the Eastern-most image 
component, along the jet axis. Slices through the core 
and along the inner jet of the aligned images are plotted 
in figure 3 and overlays of the four lowest frequencies are 
overlaid in figure 4. 

Core shifts are ignored in this approach, but these are 
expected to be small (Lobanov 1998a). Lobanov (2005, 
personal communication) predicts a core shift of 0.015 mas 
between 5 and 22 GHz for M87 assuming nominal jet 
parameters. The offset between 1.6 GHz and 43 GHz is 
0.1 mas, in line with our errors. In support of our ap- 
proach, we note the good alignment between the edges of 
the jet in the first 10 mas in the 1.6 and 4.8 GHz images. 
There is a broad component ~ 6 mas from the 1.6 GHz 
core (described in the previous section) which aligns be- 
tween the 1.6 to 22 GHz images, but it cannot be identified 
in the 43 GHz image. The 15, 22 and 43 GHz images are 
very similar and have only a small relative shift. 

It is worth noting that shifting the core site away from 
the peak of the 1.6 GHz emission does not alter the cal- 
culations of the opening angle, such as those described in 
Reid et al. 1989. An opening angle of 6.9 degrees was 
derived by Reid et al. from the increase of the jet width 
along the jet axis, and a core shift such as that described 
here therefore results in a change in the fitted core di- 
ameter at the origin from 1.7mas to 1.2mas. This core 
size matches closely the predicted size resulting from in- 
terstellar and interplanetary scattering of wlmas (Cohen 
& Cronyn 1974). 
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6.1. The spectral index 

The spectral index, a, between 1.6 and 4.8 GHz and 
4.8 to 15 GHz, after smoothing to a common resolution 
are presented in figure 5, defined in the sense S oc 
The extended emission around the 1.6 GHz image, men- 
tioned in Section 4, and not visible at higher frequencies, 
implies a spectral index of greater than 2 to the North 
and South of the region mapped. Between the higher fre- 
quencies we have a negative spectral index across most 
of the source (—0.6 ±0.4), as would expected given that 
M87 is brighter at 15 GHz than 4.8 GHz (see figure 3). 
Between 1.6 and 4.8 GHz the spectral index is more com- 
plex, with a inverted component (—1.1 ±0.4) closest to 
the core and a steep component (2.1 ±0.2) near the peak 
of the low frequency emission. This emphasizes the differ- 
ences in character of these two emission regions, and the 
analytical importance of the lower frequency, high resolu- 
tion VSOP observations. 

For the spectral fitting we used our pent achromatic 
dataset across the AGN core, and limits for the cutoff 
frequencies. Because there arc great difficulties is mea- 
suring these cutoffs we have fitted them with very large 
errors, but forcing the spectral fits to have positive curva- 
ture. We have followed the methods of Lobanov (1998b) 
to calculate the turnover frequencies as the jet develops 
along the axis. We fitted second order polynomials to the 
log(S') log{v) data, with assumed spectral cutoffs of 0.1 
and 1000 GHz, as shown in figure 6, which shows the de- 
rived maximum frequency, Vmax said the flux at that fre- 
quency, Smax- A range of cutoff frequencies were trialled 
before settling on these values. Significantly higher values 
for the upper, or smaller values for lower, cutoff distorted 
the spectrum sufficiently to discount them as possibilities. 
The two centres of emission clearly indicated in 6b, and 
suggested by the two-component image fits, place the low 
frequency peak within the coUimation zone described by 
Junor et al.(1999). We have estimated the systematic ef- 
fects of the errors in image alignment by performing multi- 
ple fits with introduced random shifts of ±0.1 mas in the 
alignment and comparing these to the plotted fits. We 
find the errors are the order of 20% for both v^ax and 

'-'max • 

Following the shocked emission approach outlined in 
Cawthorne (1991) and Lobanov (1998b) (equations 17 
and 18), wc can relate the turn-over flux density and fre- 
quency, together with the distance down the jet, to the 
magnetic field strength: 

B{r) = QA,,y-^^S;Zx 

where Co is a constant of proportionally, Vmax the turn- 
over frequency, Smax the fiux at this maximum, r the 
distance down the jet and m the power of the background 
ambient magnetic field. We use the values for the core 
together with the lower limit to the average magnetic 
field strength for the core of 22mG derived by Biretta 
et al. (1991) to determine Cq. We averaged the v^ax and 
Smax across the jet and plot these in figure 7. The fig- 
ure also plots the derived magnetic from these two values 



assuming that the ambient background magnetic field is 
proportional to r~"', where, following Lobanov (1998b) 
we consider cases with both to = 2 and to = 1. The uncer- 
tainties in Umax and Smax , and the high sensitivity of B to 
these, lead us to be cautious about the derived quantities. 
However, as a qualitative guide to the field required to 
generate the observed emission, the results are robust and 
repeatable even in our randomised datasets. Two peaks in 
the magnetic field are seen, whether ambient field varies 
as the square or linearly. The first peak occurs just be- 
fore the high frequency peak of emission and the second 
occurs just before the low frequency peak. The magnetic 
field between these two peaks is '^lOmG. This leads us 
to the speculation that core is delimited by magnetic field 
compression at the boundaries, that the field is roughly 
uniform within the core, and that the low frequency peak 
is outside this core, and in the coUimation zone. We note 
that Lobanov (1998b) derived a similar behaviour for the 
magnetic field of 3C345. 

7. Summary and Discussion 

We have presented images of M87 at 1.6 and 4.8 GHz 
from VSOP observations, and at 15, 22 and 43 GHz from 
VLBA observations three weeks later. The combination 
of lower frequency and high resolution provided by the 
space baselines allowed us to directly measure the turn- 
over frequency, rather than extrapolate it. The spectral 
index maps show distinct differences between the core and 
the component ^2 mas from the core. We have used the 
five frequencies to fit the spectrum of the source, pixel 
by pixel. From this we derive the frequency of maximum 
emission and the flux at that frequency. Wc have combine 
these values to calculate the variation of the magnetic flcld 
strength across the core and into the jet. Taking a reason- 
able set of assumptions, the magnetic field peaks at the 
outer edges of the core, and is fiat across it. 

Applying a {u,v) taper to the 1.6 GHz data and restor- 
ing with a larger beam resulted in the detection of jet 
components up to 160 mas from the core, revealing that 
these features have appeared stationary over the last 16 
years. Comparison of the 4.8 GHz image with the earli- 
est VSOP image of M87 indicates that components closer 
to the core have also appeared stationary on a 2.25 year 
timescale. 

The reported detection of TeV gamma-rays from M87 
invites comparison of these jet speeds with the sub- 
luminal and stationary components observed in other 
AGN detected at TeV energies (Edwards & Finer 2002; 
Finer & Edwards 2004). The comparison is complicated 
by the fact that M87 is an intrinsically much weaker TeV 
source (Aharonian et al. 2003). Bai & Lee (2001) consider 
the contrasting effects of M87's proximity but larger jet 
angle to the line of sight, 6. They note that if the bright 
TcV sources Mkn421 or Mkn501, both at distances of 
~150Mpc and for which they adopt 0^6°, were at the 
distance of M87 but viewed at 6* = 35°, they would be 
weak but detectable. Biretta et al. (1999) however, derive 
< 19° based on the HST observations of apparent speeds 
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of 6 c, which would increase the predicted flux (keeping 
other assumptions constant) by a factor of 20. Given all 
the uncertainties, it is at least plausible that M87 is a 
detectable (if variable; Le Bohec et al. 2004) TeV source. 

Assuming this is the case, we can then examine the 
explanations proposed for the low jet motions in TeV 
gamma-ray sources (Edwards & Finer 2002; Finer & 
Edwards 2004) for the specific case of M87. M87 shows ev- 
idence for a significant jet opening angle (Reid et al. 1989), 
which Gopal-Krishna et al. (2004) have pointed out can re- 
sult in apparent speeds much slower than the actual speed 
due to variation of boosting across the width of the jet. 
However, the effectively stationary jet motions observed 
on the mas-scale would suggest this model does not apply 
in this case. There would appear to be evidence for decel- 
erating jet model of Georganopoulos & Kazanas (2003) in 
the recent reported jet speeds close to the core at 43 GHz, 
however one must then invoke reacceleration at large core 
distances to explain the superluminal motions detected in 
HST observations. The most natural explanation for M87 
would appear to be, as discussed earlier, an invisible jet 
spine with the observed jet components being standing 
shocks in the sheath of surrounding material. 

The apparent component motions observed at 43 GHz 
thus offer the promise of studying coUimation region and, 
possibly, the transition region from significant motions to 
apparently stationary features in the parsec-scale jet. The 
planned VSOF-2 mission (Hirabayashi et al. 2004) will 
enable observations at 43 GHz with even higher resolution, 
which will enable the acceleration and coUimation of this 
important radio jet to be examined even more closely. 
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Fig. 1. Images of M87 at (a) 1.6 GHz with a resolution of 1.2 X 1.0 and contours at 5, 10, 15, 150 and 170 mjy/beam, (b) 4.8 GHz 
with a resolution of 0.7 X 0.5 and contours at 1, 2.5, 5, 10 and 100 mJy/beam, (c) 15 GHz with a resolution of 0.9 X 0.5 and contours 
at 10, 20, 40, 200 and 400 mJy/beam, (d) 22 GHz with a resolution of 0.6 X 0.3 and contours at 10, 20, 40, 200 and 400 mJy/beam, 




Fig. 2. M87 at 1.6 GHz (a) and 4.8 GHz (b) as observed by the VLBA alone and are tapered to 50% at SOMA. The 1.6GHz 
resolution is 5.9 X 2.6 and the 4.8 GHz resolution is 2.9 X 1.9. The contour levels start at 0.2 and 0.1 mjy/beam, respectively, and 
double with each level. The images have been rotated by —20.62° so that the jet runs left to right and and the 4.8-GHz image has 
been aligned so that the peak of the first emission along the jet axis coincides. Components Nl, N2 and M/L are marked. The beam 
sizes are shown in the lower right. 
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Fig. 3. Slices through the jets at five frequencies after smoothing and alignment. The contour levels start at 15 mJy/beam and 
double with each level. The four images at (b) 4.8 GHz (c) 15 GHz (d) 22 GHz and (c) 43 GHz have been shifted so that the peak 
of the first emission along the jet axis coincide with that of the 1.6 GHz (a). 



No. 



Spectral Properties and Jet Motions in M87 





Fig. 4. Overlays of the shifted images for the frequencies, (a) 1.6 GHz (greyscale and contours) and 4.8 GHz (contours) at 
1.2 X 1.0 mas, and (b) 15 GHz (greyscale and contours) and 43 GHz (white contours) at 0.9 X 0.5 and 0.3 X 0.2 mas. The con- 
tours are 5, 10 and 15 mjy/beam for 1.6 GHz; 2, 4, 8 and 16 mjy/beam for 4.8 GHz; 5, 10, 20, 40 and 80 mjy/bcam for 15 GHz; 
10, 20, 50 and 100 mjy/beam for 43 GHz. The beam sizes are shown in the lower left. 
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Fig. 5. Spectral index between (a) 1.6 and 4.8 GHz, and (b) 4.8 and 15 GHz. Images have been smoothed to 1.2 X 1.0 mas, with a 
position angle of 81°. The outline of the beam is in the bottom left of the image. Contours at every 0.5 of the spectral index are 
overlaid. Three 4.8 GHz contours also have been plotted at 2.5, 5 and 10 mjy/beam. 
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Fig. 6. a) Turnover frequency calculated from the pentachromatic data with contours at 3, 4, 5, 10, 15, and 20 GHz and b) brightness 
at that frequency with contours at 5, 10, 15, 20, 25 and 35mjy/bcam. for M87 based on fit to log(S'), log(i^). A lower cutoff of 0.1 
GHz and a high cutoff of 1000 GHz was used in the fitting. 
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Fig. 7. The one dimensional average across the jet centre of Umax, the turnover frequency in GHz, and Smax, the brightness at 
that frequency in Jy/beam (scaled up by a factor of 100) and the magnetic field in mG, assuming an underlying ambient field of 

Bamb ^ r-^, or (X r"^ 



